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In the lattice CP^" 1 models we studied the problems related to the measure of 
observables closely connected to the dynamically generated gauge field, such as the 
topological susceptibility and the string tension. We perfomed numerical simulations 
at N = 4 and N = 10. In order to test the universality, we adopted two different 



I. 

Scaling and universality tests led to the conclusion that at N = 10 the geometri 
cal approach gives a good definition of lattice topological susceptibility. On the other 
hand, N = 4 proved not to be large enough to suppress the unphysical configurations, 
called dislocations, contributing to Xt ( a ^ least up to £ ~ 30 in our lattice formula- 



tions). We obtained other determinations of Xt by the field theoretical method, wich 
relies on a local definition of the lattice topological charge density, and the cooling 
method. They gave quite consistent results, showing scaling and universality. 

The large- N expansion predicts an exponential area law behavior for sufficiently 
large Wilson loops, which implies confinement, due to the dynamical matter fields 
and absence of the screening phenomenon. We determined the string tension, without 
finding evidence of screening effects. 

PACS numbers: 11.15 Ha, 11.15 Pg, 75.10 Hk 
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I. INTRODUCTION 



Two-dimensional CP^ -1 models play an important role as a theoretical laboratory for 
testing non-perturbative analytical and numerical methods in a confining, asymptotically 
free quantum field theory. A pleasant feature of these models is the possibility of performing a 
systematic 1/N expansion around the large N saddle point solution. Indeed, most properties 
of CP^ -1 models have been obtained in the context of the 1/N expansion JE@Jj|. 

An alternative and more general non-perturbative approach is the simulation of the 
theory on the lattice. Recently there has been considerable interest in simulations of lattice 
CP*" 1 models §1,1,0,11,0,0. 

The purpose of the present paper is that of presenting rather complete numerical results 
concerning the CP 3 and the CP 9 models. We especially analyze the problems related to the 
measure of observables closely connected to the dynamical gauge field, such as the topological 
susceptibility and the string tension. 

A troublesome point in the lattice simulation technique is the study of the topological 
properties. Measuring the topological susceptibility Xt from discrete configurations proves 
to be a non-trivial task. 

It is well known that in the 2-d 0(3) a model or CP 1 model, geometrical definitions 
of topological charge are plagued by the presence of dislocations JT2] , |T3"| ], i.e. topological 
structures of the size of one lattice spacing, whose unphysical contribution to Xt does not 
vanish in the continuum limit. As a consequence, the topological susceptibility derived from 
these definitions does not show the expected scaling behavior. 

The general belief is that for higher values of N the above problems should disappear. 
But the situation for the CP 3 model appears still problematic. There have been recent 
attempts |]5|j7|j8|,p!0[| to determine Xt by using the geometrical method, but the results are not 
consistent. The authors of Ref. || and Ref. [Hj both claim to observe scaling but they find 
different values for Xt (about a factor two of difference). Since they use two different lattice 
formulations, the geometrical measure would violate universality. Furthermore, the authors 
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of Ref . |7j , using the same lattice formulation as in Ref . |L0| , do not even see scaling. 



An alternative approach relies on a definition of the topological charge density by a local 
polynomial in the lattice variables. Local operator definitions are not affected by the dislo- 
cation problem but unavoidably lead to mixing with lower and equal dimension operators 
and to the need of subtracting perturbative tails and performing finite renormalizations |TJ| . 



Neverthless, this method allowed for the determination of the topological susceptibility of 
the CP 1 or 0(3) a model 



Another important property of the CP^ -1 models is the appearance of a linear con- 
fining potential between non gauge-invariant states. The large-A expansion predicts an 
exponential area law behavior for sufficiently large Wilson loops. The persistence of the 
area law at large distance would imply the absence of the screening phenomenon due to the 
dynamical "matter fields" . The point we wish to clarify is whether at finite N the screening 
phenomenon is recovered, or the large- N prediction is confirmed. 

This paper is organized as follows. 

In Sec. [TI] the lattice actions adopted for numerical simulations are presented and the 
lattice definitions of physical observables are introduced. 

In Sec. |TJ we discuss specific features of the simulations of the CP 3 and CP 9 models, 
and present the corresponding numerical results. 

In Sec. [TV] problems related to the evaluation of the topological susceptibility are carefully 
analyzed. 

In Sec. [V] we discuss the determination of the string tension from the Wilson loops. 

II. LATTICE FORMULATION 

We choose to regularize the theory on the lattice by considering the following action: 

S s = -N(3 J2 (zn+^z An, M - 2) , (1) 

where z n is an A-component complex scalar field, constrained by the condition 
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(2a) 



and A ni u is a U(l) gauge field satisfying 



(2b) 



We also considered its tree Symanzik improved counterpart [ITS 
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gSym = _ Nf3 



^ ] (^Zn+fiZ n ^n,[j, ~\~ Z n Zn+fi^n,fj, 2^ 



n,fi 



5E( 



(3) 



Tests of rotation invariance and stability of adimensional ratios of physical quantities showed 
that the above actions lead to scaling for rather small correlation lengths |3J]. Comparison 
of measurements performed using these two actions will provide a check of universality, 
implying that the two actions are different regularizations of a unique quantum field theory. 

Since the two actions are linear with respect to each lattice variable, it is easy to construct 
efficient local algorithms based on overrelaxation procedures. In our simulations we employed 
algorithms consisting in efficient mixtures of over-heat bath [16] and microcanonical fl7[ 



algorithms. The detailed description of this simulation algorithm with a discussion of its 
dynamical features is contained in Ref. ||]. 

An important class of observables can be constructed by considering the local gauge- 
invariant composite operator 



Pij(x) = Zi(x)Zj(x) 



(4) 



and its group-invariant correlation function 



G P (x) = (TrP(x)P(0)) c 



(5) 



The standard correlation length £ w is extracted from the long-distance behavior of the 
zero space momentum correlation function ("wall- wall" correlation). The expected large- 
distance behavior, including periodic boundary condition effects, is 



G w (x) 
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(6) 



Moreover we measured the "diagonal wall- wall" correlation length obtained by summing 
on the correlations between points located on two distinct parallel lines oriented at 45° with 
respect to the coordinate axes, whose large distance behavior should be 



G A (x) * ^ 



x 
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> x » £ d 



(7) 



In practice, £ w , A w , and A d will be obtained by fitting the data for G w and G d by the 
functions (0) and (0), using all the values of x larger than a value x min to be determined. 
The comparison between £ w and £d provides a test of rotation invariance. Indeed in the 
scaling region, rotation invariance implies £ w = £d and A w = A&. Both £ w and £d should 
reproduce in the continuum limit the inverse mass gap. 

An alternative definition of the correlation length £q comes from considering the second 
moment of the correlation function Gp. In the small momentum regime we expect the 
behavior 

Z P 



G P (k) 



Co 2 + k 2 ' 



(8) 



where Gp(k) is the Fourier-transformed of Gp(x). The zero component of Gp(k) is by defi- 
nition the magnetic susceptibility Xm- On the lattice we can use the two lowest components 
of Gp(k) to obtain the following definition of £g<: 



4sin 2 7r/L 



gp(0,0) 
G P (0, 1) 



(9) 



In the scaling region the ratio £g/£w must be a constant, scale-independent number. The 
large- N expansion predicts WE 
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when iV — > 00, while for N = 2 the ratio is equal to 1 within 1% |4]]. 

The quantity Zp = Xm£c; 2 is related to the renormalization of the composite operator P. 
Its dependence on (3 can therefore be determined by renormalization group considerations. 
One finds that 



z P = c/r 2 i + o (^j , (ii) 

where c is a constant independent of the regularization scheme and therefore of the lattice 
action. In the large N limit it turns out to be 



3 

2^ 



8.5414 „ f 1 
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(12) 



N \N 2 , 

In Ref. [|J the quantity A G = Zp£ w was introduced. The adimensional ratio A w /Aq is 
another scheme independent quantity, which is approximately equal to 1 in the CP 1 (or 0(3) 
a) model H and goes to zero in the large- N limit because the zz state becomes deconfmed. 

Rotation invariance and stability of adimensional physical quantities characterize the 
scaling region. Asymptotic scaling is only needed to extract the A parameter of the lattice, 
and to check the predictions of the perturbation theory around the critical point. Of course, 
the requirement of scaling is weaker than that of asymptotic scaling, which is expected to 
be accurately testable only much closer to the critical point. On the other hand, theoreti- 
cally and for all numerical experiment purposes, the scaling property is already sufficient to 
simulate the physics of the continuum. 

Asymptotic scaling requires the ratio of any dimensional quantity to the appropriate 
power of the two-loop lattice scale 

A L oc {2tc(3) 2/n exp(-27r/3) (13) 

to go to a constant as f3~ l — > with a linear dependence on Furthermore, the ratio of 

the A parameters of two different regularizations is determined by a one loop calculation in 

perturbation theory. The ratio of the A parameters of the actions ([!]) and (0) is Q 

A Sym /O 444 \ 

I.345 exp(^iP). (14) 



A g \ N 

In the following we will also consider a modified bare coupling extracted from the energy 



j~9| , |2"0|| . The perturbative expansion of the internal energy is 

J_ 2A-1 fl_ 
2(3 + 16N 2 (3 2 + \(3 3 



B=l. + %^ + o(*) (15) 



for the action S g . Using the first term we can define 
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(16) 



/3e may be used as an alternative bare coupling. The first two terms of the perturbative 
expansion of the /5-function are universal, therefore the asymptotic scaling function in this 
new scheme is still given by Eq. ([131). Substituting (3 with (3 E in the two loop formula 
( [13]) should be equivalent to a resummation procedure which may improve the asymptotic 
behavior. From the Eq. (fL5|) we obtain the ratio of A#, the A parameter of the j3 E scheme, 



and A g : 
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The same can be done within the action S^ ym . In this case the perturbative expansion of 
the internal energy is 



1 O.lOOliV- 0.0589 1 
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(18) 



and therefore we find 

"tt(0.400A^ - 0.236) 
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(19) 



III. SIMULATIONS 

We performed Monte Carlo simulations of the CP 3 model in the two formulations (|]) 
and @ for several values of (3 corresponding to correlation lengths £ up to 30 lattice units. 
A summary of the runs is presented in Table |. There the integrated autocorrelation time 
of the magnetic susceptibility r^™ is also reported. In the following we will also show some 
data for the CP 9 model. Most of them were obtained from the simulations presented in Ref. 
Q], where the details of the runs were given. 

The data taken on lattices of different size can be used to extract the finite size scaling 
functions of the correlation length £q an d of the magnetic susceptibility for the CP 3 model. 
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In the scaling region the finite size scaling functions must be universal, that is independent 
of P and of the lattice formulation, in that they should reproduce the continuum physics in 
a periodic box. In Figs. [I] and we plot respectively fe G = £ g ,l/£g,oo and f Xm = Xm,L/Xm,oo 
versus z = L/^q for the two actions ([[]) and (0). For both actions we chose a value of (3 
corresponding to a correlation length £ ~ 4. The finite size scaling functions were obtained 
by approximating infinite lattice quantities with the corresponding values measured on the 
largest lattice available. The universality with respect to the lattice action is fully satisfied. 
From Figs. [I] and [| we also learn that z ~ 7 is a safe value where the finite size effects are 
smaller than 1%. The finite size scaling functions for the CP 9 model were shown in Ref. [|J. 

Regarding the CP 3 model, the data for the different definitions of correlation length, for 
the correlation function coefficient A w , and the ratio A w /Aq are reported in Table [TJ. The 
fits to G w and Gd were performed choosing x min ~ 2£ w ; fits using a larger x m i n gave consistent 
results. The ratios of these different definitions were analyzed by using the jackknife method. 
The models defined by S g and S^ ym enjoy rotation invariance and stability of adimensional 
physical quantities for all values of the correlation length considered. The correlation lengths 
£g an d the ratios £g/£w for the CP 9 model are reported in Table [TV[ By fitting with a 
constant the data of the above adimensional ratios, we obtained the results in Table [V|. 

£ w should reproduce in the continuum limit the inverse mass of the lowest positive parity 
state belonging to the adjoint representation. We also looked for other states, either excited 
states in the adjoint positive parity channel, or states in the other channels, the adjoint odd 
channel, the singlet even and odd channels. We did not find evidence of such states for the 
CP 3 and the CP 9 models. 

Data for the constant c of Eq. ([n]) are reported in Tables [HI] and |IV| respectively for the 



the CP 3 and the CP 9 models. Data show scaling and the two actions give very close values. 
The small discrepancies can be imputed to the non-universal terms of order f3~ l in Eq. (|Tl~D. 
We note also that the approach to scaling is slower for the action S g . 

We checked the asymptotic scaling, according to the two loop formula /(/?) = 
(27r/3) 2/Jv exp(-27r/3), by analyzing the quantity M G /A g = [£g/(/3)] _1 - 
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To begin with, in Fig. ^| we show data for the CP 1 model, which were taken by using the 
action S g ||, If this lattice CP 1 model belongs to the universality class of the 0(3) a model, 
M G /A g must tend to the asymptotic value 36.5, according to the exact result pl| . M G /A g 
appears to be constant within the errors for the largest values of (3. Its value (approximately 
47) is far from the asymptotic one. However this is not a problem, since field theory predicts 
an extremely slow approach to asymptopia for quantities like Mo/A g . In that region of 
(3, corresponding to correlation lengths from about 5 to 30 lattice spacings, the /3-function 
is well approximated by the two loop formula (indeed Mo/Ag is constant within errors of 
~ 3%) but its integral does not (the discrepancy is about 30%). It is a sort of pre-asymptotic 
region. 

The situation becomes better if we use the (3e scheme. In Fig. ^| we plot also 

M G , M g A e 

at^a^a; (20) 

where the ratio A^/A g is obtained by using Eq. fll7|) . Now data approach the correct value, 
represented in Fig. [3] by the continuous line. Is it only an accident? 

In Tables [TJ and [TV] we report data of M G /A g and M G /A g \ E for the CP 3 and the CP 9 
model. We show them respectively in Fig. [| and Fig. []. We use the A ratios given in 
Eqs. QT4] ) and (|I1]) to report all data in terms of A g . Again the (3e scheme improves the 



asymptotic scaling test. The two (3e scheme evaluations derived from the two actions S g 
and Sf ym show good agreement. As for the CP case, their value is different from those 
obtained with the standard schemes. 

IV. TOPOLOGICAL SUSCEPTIBILITY 
A. Introduction 

The topological charge density of a complex spin field z is 



q(x) = 7T- £ tw D^zD u z , (21) 

Z7T 



The topological susceptibility is denned as the correlation at zero momentum of two q(x) 
operators: 

X t= [ d 2 x(q(x)q(0)). (22) 



The large- N predictions concerning the topological susceptibility are 



^C = ^ + 0(N-^), (23) 



and 18 



^ = ^(l-^)+0(^). (24) 

Different methods have been proposed to calculate Xt on the lattice. The geometrical 
definition uses an interpolation among discrete lattice variables to assign an integer topo- 
logical charge to each lattice configuration. While for large N this definition is expected 
to reproduce the physical topological susceptibility, at low iV Xt could receive unphysical 
contributions from exceptional configurations, called dislocations, i.e. topological structures 
of the size of one lattice spacing. The dislocation contributions may either survive in the 
continuum limit, as it happens for some lattice formulations of the CP 1 or 0(3) a model, or 
push the scaling region for Xt to very large /3 values. 

Another approach relies on a definition of topological charge density by a local polynomial 
in the lattice variables. Local operator definitions are subject to mixing with lower and 
equal dimension operators and to finite renormalizations, which must be evaluated in order 
to extract Xt- 

A third method consists in measuring Xt on an ensemble of configurations cooled by 
minimizing locally the action. 

Any sensible definition of a lattice observable must show the correct scaling within each 
lattice formulation of the theory and universality among the determinations obtained with 
different lattice actions. 
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B. The geometrical definition 



The geometrical definition of the topological charge is W2 



ql = -L Im{ln[TrP n+/1+I/ P n+/1 P n ] 

+ ln[Tr P n+ „P nW P n ]}, fi^u. (25) 

Introducing the quantity # n , M = arg {z n z n+/J }, one easily obtains 

Qn = ~ A — e iiv(.@n,u, + @n+ii,v ~ On+v,u ~ @n,v) ■ (26) 
47T 

The periodic boundary conditions make the geometrical topological charge of each lattice 
configuration, Q g = J2 n 9n> integer. The topological susceptibility should then be extracted 
by measuring the following expectation value 

Xl = rr (W § ) 2 > • (27) 



V 

In Table |VI] we report the data of Xt f° r the CP 9 model. Using S g the approach to scaling 
is slow, instead for S g ym a better behavior is observed. For the leading scaling violation 
term must be 0(ln£/£ 2 ) when £ — > oo |T^]. Instead for the tree Symanzik improved actions 
the leading logarithm corrections are absent, and scale violations are 0(£~ 2 ) |23| . Assuming 
that the scaling violation term proportional to ln£/£ 2 is already dominant in our range of 
correlation lengths, we extrapolate data of Xt f° r the action S g . In Fig. || we plot Xt versus 
ln^c/^ 2 ,. A fit gives 

Xt = 0.0174(12) , 6 = 0.068(12), (28) 

where b is the coefficient of the ln^Q/^G term. The fitted value of Xt is in agreement with the 
value of Xt obtained with the action S g ym , which is Xt = 0.0176(9) at £g = 5.19(3). We then 
conclude that for the CP 9 model Xt is a good estimator of the topological susceptibility. 

The situation for the CP 3 model appears more problematic. Our data of Xt f° r the CP 3 
model are reported in Table |VT| and shown in Figs. and For both actions an apparent 
scaling is observed but data clearly violate universality. 
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C. The field theoretical method 



The field theoretical approach relies on a definition of topological charge density by a 
local polynomial in the lattice variables having the correct classical continuum limit 

q L (x) -*a 2 q(x) + 0(a 4 ) (29) 

(a being the lattice spacing). In order to determine Xt, the correlation at zero momentum 
of two q L (x) operators Xt * s calculated: 

Xt L = (E?W(0)) = ^((g L ) 2 ), (30a) 
Q L =^q L (x). (30b) 

X 

Xt is connected to Xt by a nontrivial relationship, since the presence of irrelevant operators 
of higher dimension in q L (x) induces quantum corrections. The classical continuum limit of 
q L (x) must be corrected including a renormalization constant Z(f3) pl[ . Other contributions 
originate from contact terms, i.e., from the limit x — > in Eq. (30). These contact terms 
appear as mixings with the trace of the energy-momentum tensor S(x) and with the identity 
operator J, which are the only available operators with equal or lower dimension. Therefore 
the relationship between the lattice and the continuum topological susceptibility takes the 
form 

Xt L (/5) = a 2 Z({3) 2 Xt + a 2 A({3) (S (x)) np + P((3)(I) + 0{a A ) . (31) 

np denotes the nonperturbative part (i.e., the perturbative tail must be subtracted). Z(f3), 
P{/3), and A((3) are ultraviolet effects, since they originate from the ultraviolet cutoff- 
dependent modes. They can be computed in perturbation theory as series in 

The field theoretical method consists in measuring Xt(P) by a standard Monte Carlo, 
evaluating Z((3), A((3) and P(/3), and using Eq. (|3T| ) to extract Xt- 

The requirement ( p9"D does not uniquely determine the lattice operator. Different lattice 
versions can be found and all of them should give the same physical result for Xt, instead 
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the renormalization functions Z(/3), P(/3), and A(f3) are lattice operator dependent. We 
considered two versions of lattice topological charge density operator. 



= "^E^Tr [p(x)A^P(x)A^Pi 
where is a symmetrized version of the finite derivative: 

AUp(x) = MP(x+»)-P(x-»)}. 



X, 



(32) 



(33) 



The second lattice operator is 



9ft*) = -^E^ Tr [P{x)AfP{x)^P{ 



where A^ 2 ^ is another version of finite derivative: 



(34) 



Afp{x) = |[P(x+/i) - P{x-p)\ - ±[P(x+2fi) - P(x-2n)\ . 



(35) 



g^(x) is a Symanzik tree- improved version of q\(x). 

We applied the field theoretical method to determine the topological susceptibility of the 
CP 3 model. 

In the following we will neglect the contribution of the mixing with S(x). This assumption 
is supported by a perturbative argument: the perturbative series of A(/3) starts with a f3~ 3 
term. It will be a possible source of systematic error in our calculations. 



D. The heating method 

In order to estimate the renormalization functions in Eq. (^l|) nonperturbatively, we 



applied the method proposed in Ref. p5| . We start from a configuration C carrying a 
definite topological charge Qo which is an approximate minimum of the lattice action (in 
this sense we will call it a "classical" configuration). We heat it by a local updating procedure 
in order to introduce short-ranged fluctuations, taking care to leave intact the background 
topological structure. We construct ensembles C^ ^ of many independent configurations 
obtained by heating the starting configuration Cq for the same number n of updating steps, 
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and average the topological charge over the ensembles. If £ ^> a, there should exist an 
intermediate range of n where fluctuations of length I ~ a are thermalized at the given value 
of (3 and reproduce the renormalization effects, while fluctuations at the scale I ~ £ are off 
equilibrium and still determined by the initial configuration. The average of Q L = J2 X Q L (x) 
over the configurations in this range of n should be approximately equal to Z{0) Q . 

We can also start from a constant configuration (with Q = 0) and construct other 
ensembles Cj® of configurations. We should find an intermediate region of n where the 
measure of Xt gives an estimate of the mixing P(/3) with the identity operator which, being 
a short-ranged effect (due to the fluctuations at I ~ a), is expected to be independent of the 
physical topological background structure. 

If we plot the values Q L averaged over and the values of Xt averaged over 

as functions of n, we should observe plateaus in correspondence of the above-mentioned 
intermediate ranges. The characteristics (starting point and length) of the plateaus are de- 
termined by the phenomenon of critical slowing down. The renormalization functions are 
determined by short-ranged fluctuations, which we do not expect to be critically slowed 
down; therefore the starting point of the plateaus should be independent of (3. On the other 
hand, the end point of the plateaus is reached when the Monte Carlo procedure changes 
the long-ranged modes that determine the topological properties, and critical slowing down 
should strongly affect these modes; therefore the length of the plateaus should be (3 de- 
pendent. This behavior is essential for the existence of an intermediate range of n where 
the renormalization effects can be measured: indeed the success of the present method for 
estimating Z(f3) and P(/3) strongly relies on the distinction between the fluctuations at 
distance I ~ a, contributing to the renormalizations, and those at I ~ £ determining the 
relevant topological properties. The fluctuations at I ~ a are soon thermalized, whereas the 
topological charge thermalization is much slower. 

In order to check that heating does not change the background topological structure of 
the initial configuration, after a given number of heating sweeps we cool the configurations 
(by locally minimizing the action) and verify that the cooled configurations have topological 
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charge equal to Qq. 

We used as an updating procedure a 20-hit Metropolis algorithm (tuned to 50% accep- 
tance), which gives a sufficiently mild heating. 

This method has been already applied to determine the topological susceptibility of the 
CP 1 or 0(3) a model Consistency of the direct measures of Z(j3) and P(f3) with 

the corresponding perturbative computation has been shown in Ref. fl4| within a lattice 
formulation of the 0(3) a model. 

We construct the initial configuration carrying topological charge Qq = 1 ("lattice in- 
stanton") starting from a discretization of the continuum SU(2) instanton: 

X! - X! - i(x 2 - X 2 ) 

z 1 {x) = 



p 2 + (x 1 - Xi) 2 + (x 2 - x 2 ) 2 



z 2 {x)= 



p 2 + (xi - Xi) 2 + (x 2 - x 2 ) 2 
Zi(x)— 0, i = 3, ...N , 



z(x+p)z{x) 

X ^ X)= \z(x+p)z(x)\ ■ (36) 

The parameter p controls the size of the instanton, and x is its center, which we always place 
at the lattice center: x = (L/2, L/2). Starting from the configuration (j36l), we performed a 
few cooling steps in order to smooth over the configuration at the lattice periodic boundary. 
After this procedure, we end up with a smooth configuration Cq with topological charge 
Q L w 1. The geometrical topological charge of this configuration is exactly equal to 1. 

In Fig. | we plot Qi(CW)/Q ,i and Q 2 (C^)/Q , 2 , where Q;(C«) (i = 1,2) is the lattice 
topological charge Qf = ^xlii.x) averaged over the ensemble Q , Qo,i * s the topological 
charge of the starting configuration measured by the operator Q\ . The data in Fig. |9| 
were taken at (3 = 1.12 and for the tree Symanzik improved action. We see clearly a 
plateau starting from n = 7 for both operators. For n = 12 we also cooled the sample of 
configurations finding Q\ ~ Q i after a few cooling steps. This value of n is marked by a 
dashed line in Fig. || According to the above considerations, the value of Qf at the plateau 
gives an estimate of Zi(f3). We repeated this procedure for other values of [3, and for both 
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actions (P and (||). We checked also the dependence of the measure on the size of the 
instanton p (in the range of p ~ £), and of the value of the topological charge of the initial 
configuration. The behavior of Qi(C^) /Qcm is always very similar to the case reported in 
Fig. [5] The results are presented in Table |V11J| . Configurations with topological charge two 
were constructed by allocating two instantons of size p at a distance d. 

We now proceed to the analysis of the ensembles of configurations obtained by 
heating the constant configuration defined by z{x) = (0,0,0,1) and A M (x) = 1, for 

several values of (3. In Figs. |10| and |ll] we plot the average value of Xt as a function of 
the number n of heating steps for the operator and respectively for the actions ([!]) and 
(0). For every value of (3 we observe a plateau starting from n ~ 25; the plateau is longer 
for higher values of /3, as expected. After n c heating sweeps (see Table |X|) we cooled the 
sample of configurations and found vanishing Q\ in a few cooling steps. After the plateau, 
Xt increases to reach the true equilibrium value. We identify the topological susceptibility 
measured at the plateau Xt,p with the perturbative tail at the given value of (3. Since Z(f3) 
and P{(3) have their origin in the fluctuations at I ~ a, finite size corrections are of the order 
of L~ 2 and therefore negligible on our lattice. Results are reported in Table [TX| 

The values of Z{(3) and P(f3) obtained by this procedure can be inserted in Eq. ([H]) to 
extract the physical value of the topological susceptibility. The results are summarized in 
Table |X| and shown in Figs. [7] and [8[ 

All four classes of measures (2 operators x 2 actions) show scaling within the errors. The 
small discrepancies among them should be explained by the mixing with the trace of the 
energy-momentum tensor S(x) in Eq. ([31]). Indeed the function A{(3) is operator and action 
dependent and should vary slowly in the relative small range of (3 considered. The fact that 
the discrepancies are small give further support to the initial assumption of neglecting the 
contribution of the mixing with S(x) in our calculations. From it we also get an idea of the 
systematic error of our calculations. 
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E. Heating and geometrical charge 



To clarify the origin of the failure of the geometrical method, we followed the behavior 
of the geometrical topological susceptibility Xt during the heating procedure of the flat 
configuration. In Fig. 12 we plot Xti^n^) when heating with the action S g at (5 = 1.20 and 



(3 = 1.25 (each ensembles contains 1500 configurations). Again after 40 heating sweeps 
we cooled the configurations finding vanishing topological activity. Therefore the signals 
we observe in Fig. |T2| are lattice artifacts, dislocation contributions. Notice that the value 
of xf during the heating procedure is not a negligible fraction of the corresponding results 



obtained at the statistical equilibrium (see Table |V]]). Comparing the data up to n = 40 for 
the two values of /3, we do not see evidence of critical slowing down effects, which means 
that the modes responsible for the observed signal are the short ranged (of the size of one 
lattice spacing), as dislocations are supposed to be. 

In Fig. [13] we plot data for xf (C^°^) obtained by heating at (3 = 1.12 with the action 
S"g ym . The cooling check is again performed after 40 heating sweeps. As before we observe 
an apparent topological activity after a few heating sweeps, but now the signal after 40 
sweeps is a smaller fraction of the total signal measured at the equilibrium condition. This 
should indicate that the action Sg ym is less subject to the dislocation problems. 

F. Cooling method 

We performed an independent measure of Xt using the cooling method |26] , which consists 
in measuring Xt on an ensemble of configurations cooled by locally minimizing the action 
(starting from equilibrium configurations). The idea behind the cooling method is that local 
changes should not modify the topological properties of a configuration, and its topological 
content can be extracted from the cooled configuration, where the short-ranged fluctuations 
responsible of the renormalization effects, have been eliminated. 

The cooling algorithm consists in assigning to each lattice variable z n (A„ iAt ) a new value 
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z' n {X' n „) (keeping all other variables fixed) that minimizes the action. 

To determine the topological charge of the cooled configurations, we used the operator 
which turns out to be better than (x) in estimating the topological content of a 
smooth configuration. This can be seen by comparing tha values of Qo,i and Qo,2 in Table 
Vlll|. The topological susceptibility measured on cooled configurations by Eq. (30), Xt°°\ is 



seen to gradually rise up to an extended plateau. Our averages and errors are estimated on 
the plateau measurements. 

Table [VI 1] reports also data of some measurements of Xt° o1 by cooling method for the CP 9 
model. These measures are quite consistent, especially those corresponding to the longest 
correlation lengths, with the geometrical determinations. 

For the CP 3 model the results are reported in Table [V| and plotted in Figs. [F| and §. xT° l 
shows scaling and the test of universality is satisfactory. Furthermore, they are consistent 
with the measurements obtained by the field theoretical method with the operator q\ . 



G. Conclusions 

The geometrical approach gives a good definition of lattice topological susceptibility for 
N = 10. On the other hand, we showed that N = 4 is not large enough to suppress the 
unphysical configurations contributing to Xu a ^ least for £ < 30. The other methods, field 
theoretical and cooling, give consistent measures of Xt- We finally quote for the CP 3 model 
Xt£,Q — 0-06 with an uncertainty of 10-20%. 

V. THE STRING TENSION 
A. Wilson loops and finite size effects 

In the CP^ -1 models it is possible to define the (Abelian) Wilson loop 

W(C)= J] (37) 

n,fi£C 
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The large- N expansion predicts an exponential area law behavior for sufficiently large Wilson 
loops 

W(C) ~ e -°MC)~ P P(c) for A ( C ) > ^2 ? (3g) 

where a is the Abelian string tension and p is a (renormalization-dependent) perimeter term. 
This implies also that the dynamical matter fields do not screen the linear potential at any 
distance. The large- N prediction for a is ct£q = n/N. 

Starting from the rectangular Wilson loops, we can define the Creutz ratios as 

X (l m )=ln W{l > m - l)W{l - l > m) (39) 

The double ratio takes care of renormalization effects [constant and perimeter terms in 
\nW(l,m)]. It is therefore easier to extract the string tension from x(/,m). 

In principle one can also define the Polyakov line and study the correlation of two such 
lines, thus extracting the particle-antiparticle potential. In practice the signal is so small 
and noisy that one can hardly extract a physically meaningful number. 

In order to extract the string tension from our simulations, we should understand the 
behavior of the large Abelian Wilson loop of a confining theory in a 2-d finite lattice with 
periodic boundary conditions. 

To this purpose, consider a simple 2-d model whose lattice gauge field propagator Aa is 

A,(k) = l, ^ = 4sin 2 ^. (40) 

For this model we find that a rectangular RxT Wilson loop in a L x L lattice with periodic 
boundary conditions has the following form: 



1 r d 2 k sin 2 (khR) sin 2 (\k 2 T) A 



2 J (27T)2 sin 2 sin 2 ^ 

where the integral must be evaluated by eliminating the zero mode. The result of the integral 



is 



\nW(R,T)= 1 -RT(l-^ (42) 
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As expected the propagator defined in Eq. (f40~|) gives rise to a linear confining potential with 
a string tension a — 1/2, but the finite size corrections are not small. For the Creutz ratios 
f](R) = x{Ri R) we obtain (in the following we will consider only Creutz ratios with equal 
arguments) 



v {R) = x{R,R) 



1 



1 - 



/2R-1 
V L 



(43) 



From this analysis we learn that, if there are not screening effects in the CP^ 1 models, 
for a sufficiently large R the behavior of the Creutz ratios rj(R) should be 

(2R- V 2 



7](R) ~ a 



1 - 



V L 

To compare data from different lattices it is convenient to define a rescaled Creutz ratio 



(44) 



Vr (R) = V (R) 



1 



2R- 1 
L 



a 



(45) 



A large- N prediction for the behavior of the Creutz ratios can be obtained by substituting 

in Eq. (|4l"l) the following lattice regularized version of the large- N gauge field propagator: 

, -i 



c 



1 + 



(46) 



Insertion of Eq. ( ^6|) in Eq. (^1|) allows us to define the quantity r)( N >(R) and, applying Eq. 
dD, ijW(J2). 



B. Monte Carlo results 

The gauge degrees of freedom are strongly fluctuating in the numerical simulation; there- 
fore large Wilson loops are hard to measure. The action S g allows us to define improved 
estimators for operators that are linear with respect to each X n41 variable, such as the Wilson 
loops. Improved estimators can be obtained by replacing each A n . M with its average A™p in 
the field of its neighbors: 

^ imp _ / d\ n ^X n ^ exp [2(3N Re {z n+ ^z n \ n ^)\ 
n,M / d\ n ^ exp [2(3 N Re ( 

_ Zn+ii z n Ix\^^^\^n+n z n\) (4:7) 

\z n +^ z n\ IoC2f3N\z n+ ^z n \) 
20 



where Jo and I± are modified Bessel functions. 

Another way of reducing the noise is measuring the Wilson loops on cooled configurations 
|2T|j . Few cooling steps should leave intact the longe range physical quantities, such as 
the string tension, reducing the noise coming from the short ranged modes. Cooling as 
other similar techniques, smearing and fuzzy operators, provides a sequence of approximate 
improved estimators. The Creutz ratios measured on cooled configurations as function of the 
cooling step are seen to reduce the errors and give origin to a plateau, whose length depends 
on the size of the involved Wilson loops. Then the cooling procedure starts to destroy the 
physical signal. Our averages and errors are estimated on the plateau measurements. 

To begin with we present data for the CP 9 model. In Fig. |14] we show the quantities 
r](R) £q and r] T (R) £q versus r = R/(,g- Data were obtained by using the action S g and 
at (3 = 0.8 on a 60 x 60 lattice. The improved estimators defined in Eq. (^) allow good 
measures up to r ~ 2, for larger r the signal becomes too noisy. By using cooling we found 
clear signals up to r ~ 3. As Fig. |TJ| shows, the two sets of data are in perfect agreement. 
Starting from r ~ 2 the rescaled Creutz ratios show a clear plateau which is the evidence 
of the string tension. We find <t£q = 0.25(1), to be compared with the large-iV prediction 
<Jn£,g = 0.314. This is not a surprise, in that the quantitative agreement with the 1/N 
expansion can only be reached at very large N, because of the very large coefficient in 
the effective expansion parameter 6n/N that can be easily extracted from a nonrelativistic 
Schrodinger equation analysis of the linear confining potential 

(Tat concerns the longe range predictions of the large-iV expansion. We could also test 
the short distance predictions by subtracting from r/( NS) (R) the constant a = — a and 
comparing the new curve fj^(R) with data. In Eq. (|46|) as dimensional input we use the 
value of which in the large- N limit is related to £ w by the relation (|H]). The results of 
such calculations, fj^ N ' (R) and fj^ (R) are shown in Fig. [TJ] by the continuous lines. 



In Fig. [15| we test the universality. Together with the above data we plot the rescaled 
Creutz ratios r] r (R) measured by using the action Sg ym at (3 — 0.75 on a 60 x 60 lattice. 
These last data give <t£q = 0.27(2), in agreement with the previous measure. 
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In Fig. [16| we show rj(R) and the corresponding rescaled ones r] r (R) for the CP 3 model. 
Data were taken with the action 5*g ym and at /3 — 0.95 on a 60 x 60 lattice and by using 
the cooling tecnique. Again starting from r ~ 2 the rescaled Creutz ratios shows a plateau, 
which gives a string tension <t£q = 0.31(2). 

Fig. |I7] shows the rescaled Creutz ratios coming from different simulations done with 
both actions S g and S'| yni , at several values of (3 corresponding to correlation lengths from 
about 7 to 15 lattice spacings, and on different lattices. The universality is fully satisfied. 

In conclusion, we do not see evidence of screening effects (at least up to 3£) confirming 
the qualitative picture coming from the large- N expansion. 
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FIGURES 



FIG. 1. Finite size scaling of the correlation length £g for the CP 3 model. 



FIG. 2. Finite size scaling of the magnetic susceptibility Xm for the CP 3 



model. 



FIG. 3 



Asymptotic scaling test for £g m the CP model. 



FIG. 4 



Asymptotic scaling test for £g m the CP 3 model. 



FIG. 5 



Asymptotic scaling test for £g m the CP 9 



model. 



FIG. 6. Topological susceptibility versus In £g/£g m ^ ne ^P 9 model. 

FIG. 7. Summary of the topological susceptibility determinations with the action S g and for 
the CP 3 model. 

FIG. 8. Summary of the topological susceptibility determinations with the action S^ ym and for 
the CP 3 model. 

FIG. 9. Determination of the multiplicative renormalization constants Zf for the tree Symanzik 
improved action at (5 = 1.12. 

FIG. 10. Determination of the perturbative tail P of the topological susceptibility constructed 
with the operator q\ for the action Sg. Solid lines show the plateau values. 

FIG. 11. Determination of the perturbative tail P of the topological susceptibility constructed 
with the operator for the action 5| ym . Solid lines show the plateau values. 

FIG. 12. Behavior of the geometrical topological susceptibility during the heating procedure 



using S i 
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FIG. 13. Behavior of the geometrical topological susceptibility during the heating procedure 
using Sfy™ 

FIG. 14. Creutz ratios at (3 = 0.8 with the action S g . CP 9 model. The dashed and continuous 
lines are respectively fj^ N \R) and fjr N \R). 

FIG. 15. Universal behavior of the rescaled Creutz ratios. CP 9 model. 

FIG. 16. Creutz ratios at (5 = 0.95 for the tree Symanzik improved action. CP 3 model 

FIG. 17. Universal behavior of the rescaled Creutz ratios. CP 3 model. 
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TABLES 



TABLE I. Summary of the simulation runs for the CP 3 model. Asterisk marks runs for the 
Symanzik improved action (||). We use the notation "m.,7" for a stochastic mixture of microcanon- 
ical and over-heat bath updating with relative weigth 7 (see Ref. 



p 


L 


stat 


E 


£g 




'int 


0.95 


12 


100k m,4 


0.5925(4) 


3.987(14) 


20.30(8) 


6.1(2) 


0.95 


15 


100k m,4 


0.6005(3) 


4.239(16) 


22.89(8) 


7.1(2) 


0.95 


18 


100k m,4 


0.6042(3) 


4.277(19) 


23.52(11) 


7.9(3) 


0.95 


21 


100k m,4 


0.6058(2) 


4.272(22) 


23.64(13) 


8.4(3) 


0.95 


24 


100k m,4 


0.6067(2) 


4.142(23) 


23.01(12) 


7.3(3) 


0.95 


27 


100k m,4 


0.6070(2) 


4.129(25) 


22.96(11) 


7.3(3) 


0.95 


30 


100k m,4 


0.6071(1) 


4.108(25) 


22.84(11) 


7.0(2) 


0.95 


33 


200k m,4 


0.60725(8) 


4.094(20) 


22.81(7) 


6.6(2) 


1.05 


60 


100k m,4 


0.53051(7) 


7.59(8) 


62.6(4) 


15.8(8) 


1.15 


120 


100k m,4 


0.47341(4) 


15.1(3) 


189.6(2.1) 


42(3) 


1.20a 


150 


100k m,9 


0.45018(4) 


20.9(4) 


328(4) 


64(7) 


1.206 


150 


100k m,4 


0.45011(4) 


20.2(4) 


317(4) 


70(7) 


1.25a 


210 


100k m,4 


0.42944(3) 


27.9(8) 


549(10) 


Ri 110 


1.256 


210 


100k m,9 


0.42940(3) 


28.3(7) 


552(11) 


pa 120 


0.85 * 


12 


100k m,4 


0.6529(4) 


3.936(14) 


21.54(7) 


6.0(2) 


0.85 * 


15 


100k m,4 


0.6578(3) 


4.212(17) 


24.30(10) 


7.5(3) 


0.85 * 


18 


100k ro.,4 


0.6613(3) 


4.188(20) 


24.86(12) 


8.3(3) 


0.85 * 


21 


100k m,4 


0.6631(2) 


4.160(21) 


24.82(13) 


8.2(3) 


0.85 * 


24 


100k m,4 


0.6638(2) 


4.088(22) 


24.54(12) 


8.2(3) 


0.85 * 


27 


100k m,4 


0.6640(2) 


4.025(25) 


24.21(12) 


7.6(3) 


0.85 * 


30 


100k m,4 


0.6639(2) 


4.060(27) 


24.35(12) 


7.7(3) 


0.85 * 


36 


100k m,l 


0.6641(2) 


4.044(32) 


24.28(12) 


6.9(2) 



27 



0.95 * 


60 


100k m,4 


0.57665(8) 


7.36(8) 


63.8(5) 


16.6(9) 


1.00 * 


81 


100k m,5 


0.54199(6) 


10.14(13) 


106.4(1.0) 


26(2) 


1.05 * 


120 


60k m,4 


0.51178(6) 


13.5(3) 


175(3) 


40(4) 


1.07 * 


120 


160k m,4 


0.50081(4) 


15.4(2) 


215(2) 


49(4) 


1.12 *a 


150 


160k m,4 


0.47555(3) 


21.2(4) 


360(8) 


93(9) 


1.12 *b 


150 


140k m,9 


0.47550(4) 


21.0(4) 


361(5) 


95(10) 


TABLE II. Values and ratios of different definitions of correlation length, and correlation func- 


tion coefficient for the CP 3 model. 










P 


L 


£g 








A W /A G 


0.95 


33 


4.094(20) 


4.151(21) 


0.985(4) 


0.999(5) 


0.950(13) 


1.05 


60 


7.59(8) 


7.66(14) 


0.991(10) 


1.014(15) 


0.97(4) 


1.15 


120 


15.1(3) 


15.4(4) 


0.980(16) 


0.980(20) 


0.93(5) 


1.20 


150 


20.5(3) 


20.7(5) 


0.992(13) 


0.992(17) 


0.98(4) 


1.25 


210 


28.1(5) 


28.9(1.2) 


0.973(21) 


0.98(5) 


0.91(7) 


0.85 * 


30 


4.06(3) 


4.14(4) 


0.982(6) 


0.998(8) 


0.936(19) 


0.85 * 


36 


4.04(3) 


4.12(3) 


0.982(6) 


0.984(10) 


0.935(21) 


0.95 * 


60 


7.36(9) 


7.46(10) 


0.987(11) 


0.991(14) 


0.95(3) 


1.00 * 


81 


10.14(14) 


10.30(20) 


0.984(12) 


0.98(2) 


0.95(4) 


1.05 * 


120 


13.5(4) 


13.6(5) 


0.992(21) 


0.98(3) 


0.97(7) 


1.12 * 


150 


21.1(3) 


21.4(6) 


0.983(12) 


1.01(2) 


0.95(4) 
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TABLE III. Some results for the CP 3 model. 



p 




P 2 z P 


M G /A g 




M G /A g \ E 


0.95 




1.228(7) 


39.1(2) 




26.7(1) 


1.05 




1.20(2) 


37.6(4) 




28.5(3) 


1.15 




1.10(3) 


33.9(6) 




27.6(5) 


1.20 




1.10(2) 


33.4(4) 




27.9(4) 


1.25 




1.08(3) 


32.7(6) 




27.9(5) 


0.85* 




1.067(9) 


33.4(2) 




25.3(2) 


0.95* 




1.06(2) 


32.7(4) 




26.6(3) 


1.00* 




1.04(2) 


31.7(4) 




26.5(3) 


1.05* 




1.05(3) 


31.7(8) 




27.2(6) 


1.07* 




1.04(2) 


31.3(4) 




27.0(4) 


1.12* 




1.03(2) 


30.6(4) 




26.9(3) 


TABLE IV. Some results for the CP 9 model. 


P 


£g 




P 2 z P 


M G /A g 


M G /A g \ E 


0.7 


2.35(3) 


0.971(6) 


0.905(20) 


25.73(33) 


20.57(27) 


0.75 


3.31(4) 




0.860(20) 


24.63(31) 


20.46(26) 


0.8 


4.67(3) 


0.956(6) 


0.830(8) 


23.65(14) 


20.30(13) 


0.85 


6.44(6) 


0.962(10) 


0.823(12) 


23.18(22) 


20.34(19) 


0.9 


8.83(8) 


0.973(9) 


0.816(11) 


22.88(21) 


20.42(19) 


0.7* 


3.78(2) 


0.968(6) 


0.795(7) 


22.53(11) 


19.90(9) 


0.75* 


5.19(3) 


0.974(6) 


0.787(7) 


22.18(13) 


19.94(12) 
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TABLE V. Tests of scaling by fitting 


adimensional ratio data with a constant. 


N 




S 






A W /A G 


4 




CI 

s g 


0.986(3) 




0.952(11) 


4 




cSym 
^g 


0.984(4) 




0.944(14) 


10 




s g 


0.965(4) 




0.870(13) 


10 




c-Sym 
^g 


0.971(4) 




0.900(14) 




TABLE VI. Geometric and cooled topological susceptibility for the CP 3 


model. 


P 


L 


io 4 xf 


v gt2 
Xt?G 


10 4 x cool 


,cool t2 
Xt ?G 


0.95 


33 


49.5(5) 


0.083(1) 






1.05 


60 


16.2(3) 


0.093(3) 


8.6(3) 


0.050(2) 


1.15 


120 


4.01(11) 


0.091(6) 


2.70(11) 


0.062(3) 


1.20 


150 


2.11(5) 


0.089(3) 


1.44(5) 


0.061(3) 


1.25 


210 


1.06(3) 


0.084(4) 


0.77(4) 


0.062(4) 


0.85* 


36 


40.5(6) 


0.066(1) 






0.95* 


60 


13.4(3) 


0.072(2) 


9.1(3) 


0.049(2) 


1.00* 


81 


7.0(2) 


0.072(3) 


5.2(2) 




1.05* 


120 


3.7(2) 


0.068(4) 


3.2(2) 


0.058(5) 


1.07* 


120 


2.83(7) 


0.067(2) 


2.37(8) 


0.056(3) 


1.12* 


150 


1.50(4) 


0.067(2) 


1.27(4) 


0.057(2) 
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TABLE VII. Geometric and cooled topological susceptibility for the CP 9 model. 



p 


L 


io 4 xf 


Xt^G 


10 4 x cool 


, .cool tl 
At ?G* 


0.7 


42 


50.5(1.1) 


0.0279(9) 


— 





0.75 


60 


22.6(6) 


0.0248(9) 


17.8(8) 


0.0195(10) 


0.8 


60 


9.9(3) 


0.0216(6) 


9.1(5) 


0.0197(10) 


0.85 


72 


5.1(4) 


0.0213(17) 


— 


— 


0.9 


90 


2.5(3) 


0.0198(24) 






0.7 * 


42 


12.3(6) 


0.0179(9) 






0.7 * 


60 


13.2(5) 


0.0186(7) 


12.3(4) 


0.0173(7) 


0.75* 


81 


6.5(3) 


0.0176(9) 


6.1(5) 


0.0166(13) 



TABLE VIII. Measure of the multiplicative renormalization of x L > starting from an instanton 
of size p or two instantons of size p and distance d. Qg is the geometrical charge of the initial 
configuration. The estimate of Zf" are taken by averaging the data in the range of n reported in 
the column "plateau". 



p 


L 


Qg 


P 


d 


Stat 


plateau 


Qi,o 


Z{ 


Q2,0 


7 L 
Z 2 


1.15 


48 


1 


10 




1000 


8-10 


0.9908 


0.375(7) 


0.9997 


0.420(12) 


1.15 


60 


2 


8 


12 


400 


8-10 


1.9165 


0.375(8) 


1.9912 


0.417(11) 


1.15 


60 


2 


8 


16 


400 


8-10 


1.9517 


0.377(7) 


1.9971 


0.418(10) 


1.20 


48 


1 


10 




2000 


7-10 


0.9908 


0.413(5) 


0.9997 


0.461(7) 


1.20 


60 


2 


8 


20 


400 


7-10 


1.9590 


0.413(7) 


1.9983 


0.459(10) 


1.20 


48 


1 


10 




1000 


7-10 


0.9908 


0.432(6) 


0.9997 


0.481(10) 


1.07* 


48 


1 


10 




1000 


7-10 


0.9905 


0.442(6) 


0.9997 


0.494(9) 


1.12* 


48 


1 


10 




1000 


7-12 


0.9905 


0.464(5) 


0.9997 


0.518(8) 


1.17* 


48 


1 


10 




500 


7-12 


0.9905 


0.487(7) 


0.9997 


0.541(11) 
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TABLE IX. Measure of the perturbative tail Pi of xt- Data were taken on 36 x 36 lattice. 






Stat 




plateau 


n c 


10 4 Pi 


10 4 P 2 


1.15 


5000 




23-30 


28 


0.293(6) 


0.714(14) 


1.20 


6000 




24-40 


35 


0.244(4) 


0.592(12) 


1.25 


6500 




28-40 


40 


0.202(3) 


0.491(9) 


l.oU 


1UUU 




on /in 
ZZ— 4U 


/in 
4U 


U.l I \)y ( ) 


\)AZ\Z ) 


1.07 * 


3000 




24-33 


28 


0.215(6) 


0.508(13) 


1.12 * 


3000 




25-38 


40 


0.176(6) 


0.414(10) 


1.17 * 


1000 




22-40 


40 


0.140(6) 


0.33(2) 


TABLE X. Measure of xt by the field theoretical method. 




io 4 xf 


10 4 Xt 


At ?G 




i n4 ft, 2 


v ft 'V 2 

At SG 


1.15 


0.647(18) 


2.52(17) 


0.058(4) 


1.265(36) 


3.1(3) 


0.071(7) 


1.20 


0.485(12) 


1.41(8) 


0.060(4) 


0.952(22) 


1.68(14) 


0.071(5) 


1.25 


0.345(8) 


0.76(5) 


0.060(5) 


0.706(15) 


0.92(9) 


0.072(8) 


1.07* 


0.646(16) 


2.21(11) 


0.052(4) 


1.147(28) 


2.62(16) 


0.062(4) 


1.12* 


0.442(10) 


1.24(6) 


0.055(3) 


0.800(16) 


1.44(9) 


0.064(4) 
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